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Abstract. We show convergence of solutions to equilibria for quasilinear 
parabolic evolution equations in situations where the set of equilibria is non- 
discrete, but forms a finite-dimensional C 1 -manifold which is normally hyper- 
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functional as in the Lojasiewicz-Simon approach, but are of local nature. 
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1. Introduction 

The principle of linearized stability is a well-known and powerful tool for proving 
stability or instability of equilibria of nonlinear evolution equations. It is known 
to be true for large classes of nonlinear evolution equations, even for such which 
are nonlocal. The literature on this subject is large. Since here we are mainly 
interested in quasilinear parabolic problems, we only refer to the monograph by 
Lunardi [21] , and to [H . 

In this paper we will consider the following situation: suppose that for a nonlin- 
ear evolution equation we have a ^-manifold of equilibria £ such that at a point 
it* G £, the kernel N(Aq) of the linearization Aq is isomorphic to the tangent space 
of ^ at «», the eigenvalue of A is semi-simple, and the remaining spectral part 
of the linearization A is stable. Then solutions starting nearby exist globally 
and converge to some point on £. This result is well-known to specialists in the 
area of dynamical systems (where it is considered a folk theorem), but might be 
less familiar to people in the PDE community. 

The situation described above occurs frequently in applications. We call it the 
generalized principle of linearized stability, and the equilibrium u» is then termed 
normally stable. 

A typical example for this situation to occur is the case where the equations 
under consideration involve symmetries, i.e. are invariant under the action of a 
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Lie-group Q. If then u* is an equilibrium, the manifold £ includes the action of Q 
on u* and the manifold Qu* is a subset of £ . 

A standard method to handle situations as described above is to refer to center 
manifold theory. In fact in that situation the center manifold of the problem 
in question will be unique, and it coincides with £ near u*. Thus the so-called 
shadowing lemma in center manifold theory implies the result. Center manifolds 
are well-studied objects in the theory of nonlinear evolution equations. For the 
parabolic case we refer to the monographs [HI [2T] , and to the publications [7J 

nana EDiEaEgisg]. 

However, the theory of center manifolds is a technically difficult matter. It 
usually involves higher regularity of the involved nonlinearities - in particular con- 
cerning the shadowing property. Therefore it seems desirable to have a simpler, 
direct approach to the generalized principle of linearized stability which avoids the 
technicalities of center manifold theory. 

The purpose of this paper is to present such an approach. It turns out that 
the effort is only slightly larger than that for the proof of the standard linearized 
stability result - which is simple. We emphasize that our approach requires only 
C^-regularity for the nonlinearities. By several examples we will illustrate that 
our result is applicable to a variety of interesting problems in different areas of 
applied analysis. It is our belief that the approach devised in this manuscript will 
be fruitful for the stability analysis of equilibria for parabolic evolution equations 
that involve symmetries in the way described above. 

Here we would also like to mention the work in [9], where the action of a Lie 
group has been used for the stability analysis of equilibrium solutions. However, 
the approach given here is considerably more general and flexible. 

In Section 2 we formulate and prove our main result for abstract autonomous 
quasilinear parabolic problems. Theorem 12.11 implies . for instance, the main result 
in [15] on convergence of solutions for the Mullins-Sekerka problem. We also show 
by means of examples that the conditions of Theorem 12.11 are necessary in order 
to have convergence to a single equilibrium. 

In Section 3, we consider quasilinear parabolic systems with nonlinear boundary 
conditions and we show that our techniques can also be applied to this situation. 
Sections 4 and 5 illustrate the scope of our main result, as we show convergence 
towards equilibria for the Mullins-Sekerka model, and stability of travelling waves 
for a quasilinear parabolic equation. 

In Section 6 we consider the so-called normally hyperbolic case, where the re- 
maining part of the spectrum of A$ also contains an unstable part away from the 
imaginary axis. In this situation, one cannot expect convergence of all solutions 
starting near it*, but only for those initial values which are on the stable manifold. 

To cover the quasilinear case our approach makes use of maximal L p -regularity 
in an essential way. As general references for this theory we refer to the recent 
publications QUE], to the survey article [24], and also to [2l l3l l4l l8l |2T] . 

In a forthcoming paper these results are extended to the case where the bound- 
ary conditions are of relaxation type, i.e. are coupled with an evolution equation on 
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the boundary, as in [13] . Problems of the last kind are important e.g. for the Ste- 
fan problem with surface tension, see [HJ [26] , and for the two-phase Navier-Stokes 
problem with a free boundary. 

Finally, we should like to point out that the generalized principle of linearized 
stability described in the current paper can also be adapted and applied to fully 
nonlinear parabolic equations, see [27j . 

2. Convergence for abstract quasilinear problems 

Let Xq and X\ be two Banach spaces such that X\ Xq, i.e. X\ is contin- 
uously and densely embedded in Xq. In this section we consider the autonomous 
quasilinear problem 

ii(t)+A(u(t))u(t)=F(u(t)), t>0, u{0) = u . (2.1) 

For 1 < p < oo we introduce the real interpolation space X y := (Xq, Xi) 1 _ 1 / p p 
and we assume that there is an open set V C X 7 such that 

(A,F)eC 1 (V,B(X 1 ,Xo)xX ). (2.2) 

Here B(Xi,Xq) denotes the space of all bounded linear operators from X\ into 
Xq. In the sequel we use the notation | ■ \j to denote the norm in the respective 
spaces Xj for j = 0, 1,7. Moreover, for any normed space X, Bx(u,r) denotes 
the open ball in X with radius r > around u £ X. 

Let £ dV C\Xi denote the set of equilibrium solutions of (|2.1[) , which means that 

it £ £ if and only if u £ V n Xi, A(u)u = F(u). 

Given an element u* £ £, we assume that it* is contained in an m-dimensional 
manifold of equilibria. This means that there is an open subset U C M. m , £ U, 
and a C 1 -function ^ : U — * X\, such that 

• *([/) C £ and *(0) = tt*, 

• the rank of ^'(O) equals to, and (2-3) 

. A(*(C))*(C) = f(*(0)i Cet/. 

We assume further that near it* there are no other equilibria than those given by 

Le - £ dB Xl (u*,ri) = *(£/), for some r x > 0. 
We suppose that the operator A(u<,) has the property of maximal L p -regularity. 
Introducing the deviation v = u — from the equilibrium u*, the equation for v 
then reads as 

ii{t)+A v(t) = G(v{t)), OO, v(0)=v , (2.4) 
where vq = uq — it* and 

A Q v = A(u*)v + - F'(u*)v for v £ X x . (2.5) 

The function G can be written as G(v) = G\(v) + 62(1*, 1 1 ), where 

Gt(v) = (F(«.+«)-F(u,)-f (u,)w)-(i4(u*+«)-A(«.)-^'(«*)«K. 
G 2 (v,w) = -(A(ti* +v)—A(u*))w, w £ Xx, v £ V*, 
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where V* := V — u*. It follows from (|2~2j) that Gi G G 1 ^*, A ) and also that 
G 2 6 C^iy* x Xi,Ao). Moreover, we have 

Gi(0) = G 2 (0,0) = 0, G' 1 (0) = G^(0,0)=0, (2.6) 

where G' 2 and G' 2 denote the Frechet derivatives of G\ and G 2 , respectively. 

Setting ip(() = ^(C) — w* results in the following equilibrium equation for problem 

mm 

A 01 p(0 = G(V»(C)), for all C € CA (2.7) 
Taking the derivative with respect to £ and using the fact that G'(0) = we 
conclude that Aotp'(0) = and this implies that 

T u ,{£) CN(A ), (2.8) 

where T Ut (£) denotes the tangent space of £ at u*. 

After these preparations we can state the following result on convergence of solu- 
tions starting near it*. 

Theorem 2.1. Let 1 < p < oo. Suppose u* £ V D X\ is an equilibrium of i2.1\) , 
and suppose that the functions (A, F) satisfy (|2.2p . Suppose further that A(u*) 
has the property of maximal L p -regularity. Let Aq, defined in (|2.5p . denote the 
linearization of h2.1\) at u*. Suppose that is normally stable, i.e. assume that 

(i) near u* the set of equilibria £ is a C 1 -manifold in X\ of dimension m G N, 

(ii) the tangent space for £ at u* is given by N(Aq), 

(iii) is a semi-simple eigenvalue of Aq, i.e. N(Aq) © R(Aq) = Xq. 

(iv) o-(Aq) \ {0} C C+ = {z e C : Rcz > 0}. 

Then is stable in X^, and there exists S > such that the unique solution u(t) 
of 112.1)) with initial value Uq S X 1 satisfying \uq — u*| 7 < 8 exists on K_|_ and 
converges at an exponential rate in X 7 to some £ £ as t — > oo. 

Proof, (a) Note first that assumption (iii) implies that is an isolated spectral 
point of o~(Aq), the spectrum of Aq. According to assumption (iv) o~(Aq) admits 
a decomposition into two disjoint nontrivial parts with 

a(A ) = {0}Uer s , <r s cC + = {zeC: Rez>0}. 

The spectral set a c :— {0} corresponds to the center part, and a s to the stable 
part of the analytic Go-semigroup e~ Aot , or equivalently of the Cauchy problem 

W + AqW = f. 

In the following, we let P , I E {c, s}, denote the spectral projections according 
to the spectral sets a c — {0} and a s , and we set X!j := P l Xj for / 6 {c, s} and 
j G {0, 1,7}. The spaces Aj are equipped with the norms | • \j for j = 0, 1, 7. We 
have the topological direct decomposition 

Ai = AJ © Af, A = A c ffiA s , 

and this decomposition reduces Aq into ^4 = A c ® A s , where A\ is the part of 
Aq in Aq for I G {c, s}. Since a c = {0} is compact it follows that Aq C X\. 
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Therefore, Xfi and Xf coincide as vector spaces. In the following, we will just 
write X c = (X c , \ ■ \j) for either of the spaces X§ and Xf. We note that N(Aq), 
the kernel of Aq 7 is contained in X c . The operator A s inherits the property of L p - 
maximal regularity from Aq. Since cr(A s ) — a s C C+ we obtain that the Cauchy 
problem 

w + A s w = /, w(0) = 0, (2.9) 
also enjoys the property of maximal regularity, even on the interval J = (0, oo). 
In fact the following estimates are true. For any a G (0, oo] let 

Eo(o)=i p ((0,o);JC ), Ei (a) = ^((0,a);X o ) n L p ((0,a); (2-10) 

The natural norms in Ej(a) will be denoted by || ■ ||E-(a) for j = 0,1. Then the 
Cauchy problem (|2.9p has for each / G L p ((0,a); Xq) a unique solution 

W ei7 p 1 ((0,a);X s )nL p ((0,a);X 1 s ), 

and there exists a constant M such that |w| El ( a ) < Mo|/[e ( ) f° r every a > 0, 
and every function / G L p ((0, a); Xq). In fact, since a(A s — oj) is still contained in 
C+ for oj small enough, we see that the operator A s — ui enjoys the same properties 
as A s . Therefore, every solution of the Cauchy problem (|2.9p satisfies the estimate 

|e fft u>|| El( a) <Mo||e CTt /|E (a), <re[0,w], a > Q > ( 2 - n ) 
for / G L p ((0,a);X o s ), where M = M (w) for to > fixed; cf. [H Sec. 6]. 
Furthermore, there exists a constant Mi > such that 

\\e u ' t e- A > t P 3 u\\ Ma} + sup |e w *e- A «*P s ti| 7 <Mi|P s u| 7 (2.12) 

te[0,a) 

for every u G X 7 and a G (0, oo]. For future use we note that 

sup |w(t)| 7 < CofltulEjfo) for all w G Ei (a) with w(0) = (2.13) 

te[0,a) 

with a constant Co that is independent of a G (0, oo], see for instance the proof of 
[25l Proposition 6.2]. We remind that iV(^4o) is contained in X c . 

(b) It follows from the considerations above and assumptions (i)-(iii) that in fact 

N(A ) = X c and dim(X c ) = m. 

As X c has finite dimension, the norms | • \j for j — 0, 1,7 are equivalent, and we 
equip X c with one of these equivalent norms, say with | • | . Let us now consider 
the mapping 

g:UcR m ^X c , g(() := F C ^(C), ( € U. 
It follows from our assumptions that g'(0) = P c ')/''(0) : M m — > X c is an iso- 
morphism (between the finite dimensional spaces R m and X c ). By the inverse 
function theorem, g is a C^-diffeomorphism of a neighborhood of in R m into 
a neighborhood, say Bx<={0,po), of in X c . Let g^ 1 : £?x<=(0,po) — > U be the 
inverse mapping. Then g^ 1 : Bx^(0,po) — > J7 is C 1 and <7 _1 (0) = 0. Next we set 
:= ■0(.9 _1 (-t)) for a; G Px=(0,po) and we note that 

3>£C\B X c{0,p ),Xl), $(0) = 0, {$( I ) + U ,:^Bx c (0, Po )}-^nW, 



6 



J. PRUSS, G. SIMONETT, AND R. ZACHER 



where W is an appropriate neighborhood of it* in X\. One readily verifies that 

P c <S>(x) = {{P c o^)og- l ){x) = (gog- 1 )(x)=x, xGB x *(0,p ), 

and this yields $(x) = P c $(x) + P s $(x) = x + P s <f>(x) for x G B X c(0,p o ). Setting 
(j)(x) := P s $(x) we conclude that 

0eC rl (Bjc.(Q,A,),*?) J 0(0) = 0'(O) = 0, (2.14) 

and that 

{x + <Kx) + u* : a; € B x <= (0, p )} = £ H W, 
where is a neighborhood of it* in Xi . This shows that the manifold £ can be 
represented as the (translated) graph of the function in a neighborhood of it* . 
Moreover, the tangent space of £ at u* coincides with N(Aq) — X c . By applying 
the projections P l , I £ {c, s}, to equation (|2.7p and using that x+(f>(x) = ip(g^ 1 (x)) 
for x 6 -Bx c (0, po)j and that A c = 0, we obtain the following equivalent system of 
equations for the equilibria of Ij2.4|) 

F c G(x + 0(a)) = 0, P s G(x + <p{x)) = A s <f>(x), xeB Xc {0, Po ). (2-15) 
Finally, let us also agree that pq has already been chosen small enough so that 

W(x)\B(x;xfi < 1, l^)|i < N, x G B X c(0,po)- (2.16) 
This can always be achieved, thanks to (12.141) . 
(c) Introducing the new variables 

x = P c v = P c (u-u*), 

y = P s v - (j)(P c v) = P s {u - u.) - <t){P c {u - i/*)) 
we then obtain the following system of evolution equations in X c x Xq 

| x=T(x,y), x(0)=x , ^ ^ 

+ = R(x,y), y(0) = y , 

with xq — P c vq and j/q = P Sv q ~ <j>(P c vo)> where the functions T and R are given 

by 

T(x,y) = P c G(x + <f>(x)+y) 7 

R{x, y) = P s G(x + 4>{x) + y) - A a cf>{x) - <j>'{x)T{x, y). 

Using the equilibrium equations (|2.15j) , the expressions for R and T can be rewrit- 
ten as 

T(x, y) = P c (G(x + <j>{x) + y) - G(x + cj>(x))) , 

(2.1o) 

R(x, y) = P s (G(x + 4>(x) +y)- G(x + cf>(x))) - </>'(x)T(x, y). 

Although the term P c G(x + <f>(x)) in T is zero, see (|2.15j) . we include it here for 
reasons of symmetry, and for justifying the estimates for T below. Equation (|2.18j) 
immediately yields 

T(x, 0) = R(x, 0) = for all x G B X c(Q, Po), 
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showing that the equilibrium set £ of ()2.1|) near has been reduced to the set 
B X c{O lPo ) x {0} C X c x Xf. 

Observe also that there is a unique correspondence between the solutions of (|2.1|) 
close to u* in X 1 and those of (|2.17j) close to 0. We call system (|2.17j) the normal 
form of (|2.ip near its normally stable equilibrium u* . 

(d) From the representation of G and (|2.6[) we obtain the following estimates for 
G\ and G 2 : for given r\ > we may choose r = r(^) > small enough such that 

|Gi(t>i) - Gi{v 2 )\o < r)\vi - v 2 \ 1 , vx,V2 G Bx 7 (0,r). 

Moreover, there is a constant L > such that 

\G 2 (vi,w) - G 2 (v 2 ,w)\ < L\w\i \vi - v 2 \-y, w G X\, vi,v 2 G Bx^(0,r), 

\G 2 (v, wi) - G 2 {v, w 2 )\ < L r \wi - w 2 \i, Wi,w 2 G X\, v e B Xl (0,r). 

We remark that L does not depend on r G (0, ro] with ro appropriately chosen. 
Combining these estimates we have 

\G{vi) - G(v 2 )\ < (r) + L\v 2 \i)\vi - u 2 | 7 + Lr\v\ - v 2 \\ 

< C (r) + r + |u 2 |i)K -V2I1 

for all Wi,f2 G -6^(0, r) fl Xi, where Co is independent of r G (0,ro]. 

In the following, we will always assume that r 6 (0,ro] and ro < 3po- Taking 
v x = x + <j){x) + y and v 2 = x + <j){x) in (f27T9|) we infer from ([2~T6]) and ([2~T8]) that 

|i?(x,y)|o<C 1 (7 ? + r+|x + </)( : r))| 1 )|y| 1 </3|y| 1 , (2.20) 

for all x £ Bx^O,^), ?/ G 5x^(0, p)nXi and all p G (0,r/3), where (3 = C 2 {r} + r), 
and where C\ and C2 are uniform constants. Suppose that r\ and, accordingly, r 
were already chosen small enough so that 

M /3 = M C 2 (f] + r)< 1/2. (2.21) 

(e) By [24l Theorem 3.1], problem (12. 4| admits for each vq G 5j 7 (0, r) a unique 
local strong solution 

u€E 1 (o)nC([0,o];X r ) (2.22) 

for some number a > 0. This solution can be extended to a maximal interval of 
existence [0, i*). If is finite, then either v(t) leaves the ball Bx^ (0, r) at time i*, 
or the limit lim^^ v(t) does not exist in X 7 . We show that this cannot happen 
for initial values vo G i?x T (0, 5), with 8 < r to be chosen later. 

Suppose that xq G £?x c (0,^(5) and ?/o G £?x s (0, iV<5) are given, where the number 
8 will be determined later and N :— \\P c \\t3(x ) + \\P S \b(x^)- Let denote the exis- 
tence time for the solution (x(t), y(t)) of system (I2.17[) with initial values (xo, yo), 
or equivalently, for the solution v(t) of (12. 4|) with initial value vq = XQ + (f>(xo) + yo. 
Let p be fixed so that the estimates in (|2.20| hold. Set 

h := h(x ,yo) ■= sup{£ G (0,4) : \x(t)\, |y(r)| 7 < p, r G [0,t]} 
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and suppose that ti < t*. Due to (12.11|) - (12.12p and (|2.20[) we obtain 

|e wt J/|| El ( tl ) < + Mo\\e^R(x,y)\\ Eo{tl) 

< Mx\y \ 7 + M Q p\\e ut y\\ El{tl) . 

This yields with (|2~2l"l) 

le^tflEiCti) < 2Mi|y | 7 , a G [0,w]. (2.23) 

Using this estimate as well as f|2.12|) ()2.13[) we further have for t G [0, t\) 

\e at y(t)\ 7 < \e^y{t) - e^e^y^ + ^e^y^ 

<c |e <rt y(t)-e <rt e- A ' t W) | El(tl) +M 1 | I «,| 7 

< (3coAfi+Afi)|ito| 7 , 

which yields with M 2 = (3co + l)Mi, 

\y(t)\ y < Af 2 e- <rt | W )| 7 , * £ [0,*i), a G [0,w]. (2.24) 

We deduce from the equation for x, the estimate for T in (|2.20[) . and Holder's 
inequality that 

\x(t)\ < \xa\+ f \T{x(s),y(s))\ds 
Jo 

< \xo\+P I \y{s)\ids 



< 



M + /?( J~ e-^' ds) 1/P ' \\e^y\\ El{tl) 



= \x \ + pcxie^yiE^) < \x \ + M 3 \y \^ ie [0,*i), 

where M3 = M\C\/Mq and c\ = (\/[u>p']) 1 / p . Summarizing, we have shown that 
\x{t)\ + |j/(i)| 7 < |x | + (M 2 + M 3 )\y \-y for all t G [0,*i). By continuity and the 
assumption t\ < t* this inequality also holds for t = t\. Hence 

|x(ti)| + |y(ti)| 7 < \x Q \ + (M 2 + M 3 )\y \-y < (1 + M 2 + M 3 )NS < p/2, 

provided S < p/[2N(l + M 2 + M3)]. This contradicts the definition of t\ and we 
conclude that t\ =t*. 

In the following, we assume that 5 < p/[2N(l + M 2 + M3)]. Then the estimates 
derived above and 12 . 16[) yield the uniform bounds 

Hk(a)+ sup Kt)| 7 <M, (2.25) 
te[o,o) 

for every initial value G -Bx T (0, S) and every a < f*. It follows from Corollary 3.2 
in [21] that the solution v(t) of (|2.4[) exists on K + . 

(f) By repeating the above estimates on the interval (0, 00) we obtain the estimates 
\x(t)\<\x \+M s \y \^, \y(t)\ 7 < A'he-^lyoU, te[0,oo), (2.26) 



CONVERGENCE OF SOLUTIONS TO EQUILIBRIA 



for all x £ B X c{0, NS) and y £ B x -(0, NS). Moreover, 

lim x(t) = x + / T(x(s),y(s))ds =: Xoo 

t~,oo J Q 

exists since the integral is absolutely convergent. Next observe that we in fact 
obtain exponential convergence of x(t) towards Xoc, as 



\x(t) 



<0 



T(x(s),y{s))ds 



\y(s)\ids 



ds 



i/p' 



l|e" s 2/| El (oo) 



This yields existence of 

Voo ■ 



lim v(t) 

t — >oo 



<M 4 e- w *|y | 7 , t>0. 



lim x(t) + 4>{x{t)) + y(t) = x c 

t — >oo 



(^oo ) ■ 



Clearly Voo is an equilibrium for equation (|2.4| , and + u» £ £ is an equilibrium 
for (|2.2ll . Due to (|2.16|) . (|2.26|) and the exponential estimate for \x(t) — x^j we get 

\v{t) - Woc| 7 = \x(t) + <f>(x(t)) + y{t) - Woo| 7 

< |x(t) - a;oo| 7 + \<p(x(t)) - (j>(xoo)\^ + |t/(t)| 7 

<(CM A + M 2 )e- ut \y U (2 ' 27) 

< Me- ut \P s v Q - <t>(P c v )\^ , 

thereby completing the proof of the second part of Theorem 12.11 Concerning 
stability, note that given r > small enough we may choose < S < r such that 
the solution starting in Bx„(u*,5) exists on R + and stays within Bx (ii*,r). □ 

Remarks 2.2. (a) Theorem 12.11 shows . given that situation, that near it* the set 
of equilibria constitutes the (unique) center manifold for (|2.1[) . 

(b) It is worthwhile to point out a slightly different way to obtain the function 
<p used in the proof of Theorem 12. II Applying the projections P s and P c to the 
equilibrium equation (|2.7p yields the following equivalent system of equations near 
v = 

A s z = P s G(x + z), A c x = P c G{x + z), (2.28) 

with z = P S V(C) and x = P c ip((). Since G(0) = G'(0) = and A s is invertible, 
by the implicit function theorem we may solve the first equation for z in terms of 
x, i.e. there is a C 1 -function : Bx<=(0, po) — > X* such that 

0(0) = and A s (j){x) = P s G(x + <j)(x)), x£B x -(0,pa)- 

As x + <j){x) is the unique solution of the first equation in (12.28|) we additionally 
have A c x = P c G(x + 4>{x)), as well as P s ip(() = (f>(P c tp(Q) for all £ £ U. Since 
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G'(0) = we obtain A a <j/(0) = P s G'(0) = and this implies 0'(O) = 0. This 
shows that £ C A4 with A4 = {x + <f>{x) + u* : i£ (0, po)} in a neighborhood 
of it* in X[ . 

A( is a C 1 -manifold of dimension £ := dim (A c ) with tangent space T u „ (Af) = 
X c and £ is a submanifold in At. In general, £ has lower dimension than M.. Our 
assumptions in Theorem l2.1l do in fact exactly amount to asserting that £ and M. 
are of equal dimension. Since £ C M. we can then conclude that they coincide in 
a neighborhood of it*. 

(c) An inspection of the argument given above shows that in fact all equilibria 
of equation (|2.ip that are close to the equilibrium it* are contained in a manifold 
M = {x + <f>(x) + it* : x G Bx-(0, po)} such that 0(0) = $6'(0) = 0, with no 
additional assumptions on the structure of the equilibria. To see this, let us once 
more consider the equation 

A s z = P s G{x + z), x£X c , zeXf. (2.29) 

Clearly, x = z = is a solution. Exactly as in the remark above, we can solve 
(I2.29P by the implicit function theorem for z in terms of x, obtaining a C 1 -function 
4> : B X c(0,po) —> Xf with 0(0) = cf>'(0) = 0. If v G X% is an equilibrium for the 
evolution equation (|2.4j) close to 0, then the pair x — P c v, z = P s v necessarily 
satisfies equation (|2.29j) . and therefore lies on the graph of 4>. 

(d) We illustrate by means of examples that convergence to equilibria fails if 
one of the conditions (i)-(iii) in Theorem 1 2 . 1 1 does not hold. 

Example 1. Consider in G := R 2 \ {0} the ODE system 

x=(x + y)(l- y /x 2 +i/ 2 ~ ), , 2 3Q > 

y=(y- x)(l - ^Jx 1 + y 2 ). 

In polar coordinates (|2 . 30[) reads as 

r = —r(r — 1), 
9 = r - 1, 

thus the set of equilibria £ of (|2.30p in G is the unit circle, and for any initial value 
(^OjZ/o) G G we have r(t) — > 1 as t — > oo. Since the phase portrait is rotationally 
invariant we may restrict the stability analysis for £ to one equilibrium, say it* = 
(0, 1). Denoting the right-hand side of ([230)1 by F(x, y), we have F G C 1 (G), and 

A = -F'(u*) = 

The eigenvalues of Ao are and 1 with eigenvectors (1, 0) and (1, 1), respectively. 
Thus is semi-simple, and N(Aq) coincides with the tangent space T Ut (£). Con- 
sequently, is normally stable, and hence we can apply Theorem 1 2. II to conclude 
that each trajectory converges to some point on the unit circle as t — > oo. It is 
readily seen that the trajectories satisfy the relation 9(r) = cq — lnr for some 
appropriate constant cq, and this confirms that 9(r) converges as r — > 1. 



1 
1 
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Example 2. In this example, we consider in G := M. 2 \ {0} the ODE system 

x = -xj ^F+W - If - yj ^+y 2 ~ l) m , . . 

V = -y{\/x 2 + y 2 - l) 3 + x( y/x 2 + y 2 - l) m , [ ' ' 

with m = 1. In polar coordinates (|2 . 31[) reads as 

r = — r(r — 1) , 
= ( r _!)"». 

Again, the set of equilibria £ of (|2.31[) in G is the unit circle. As above, we may 
restrict the stability analysis for £ to one equilibrium, say = (0, 1). Denoting 
the right side of (|2.31[) by F(x, y) we obtain (in case m = 1) 



A = -F'{u* 



1 




Clearly, {0} is an eigenvalue of Aq with algebraic multiplicity 2, and N(Aq) = 
span{(l,0)}. Therefore, the eigenvalue {0} has geometric multiplicity 1 and alge- 
braic multiplicity 2. So we have the following situation: 

(i) £ — {(x,y) : x 2 + y 2 = 1} is a smooth manifold of dimension 1 in M 2 , 

(ii) the tangent space of £ at it* is given by N(Aq), 
(hi) {0} is not semi-simple, 

and hence condition (iii) of Theorem 12.11 is not satisfied. We will show that the 
trajectories of system (|2.31l) still converge towards the unit circle, but will spiral 
around the circle at increasing speed as r — > 1. This can be seen as follows. First 
we observe that V(x, y) := (r — l) 2 with r = ^/ x 2 + y 2 is a Lyapunov function for 
system (|2.31l) . since for every solution (x,y) of (|2.31[) we have 

^-V(x,y) = 2r(r - 1) = -2r(r - l) 4 < 0. 
at 

So r(t) — ► 1 as t — ► oo for every solution. On the other hand one verifies that the 
trajectories satisfy the relation 9{r) = c + ln(|r — l|/r) + l/(r — 1). This shows 
that all trajectories spiral around £ with increasing speed, in clockwise direction 
as r y 1, and in counter-clockwise direction as r \ 1. 

Example 3: Here we consider system (|2.31|) with m = 2. This example is similar 
to the one in [51 p. 4] . In this case we have 



A = -F'(u*) = 








Clearly, {0} is now an eigenvalue with geometric multiplicity 2, and N(Ao) = M 2 . 
So condition (ii) of Theorem [2J] is not satisfied. The function V from the previous 
example is again a Lyapunov function, and this yields r(t) — > 1 as t — > oo. The 
trajectories satisfy 9(r) = cq — ln(|r — l|/r), showing that they spiral counter- 
clockwise with increasing speed around the unit circle as r — * 1. 
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3. QUASILINEAR PARABOLIC PROBLEMS WITH NONLINEAR BOUNDARY 

CONDITIONS 

The analysis in the previous section applies in particular to quasilinear para- 
bolic systems of partial differential equations with linear autonomous boundary 
conditions. In this section we show how this can be extended to the case where 
also the boundary conditions are nonlinear. For this purpose, let ft C M™ be an 
open bounded domain with boundary dft G C . The outer normal at a point 
x E <9ft will be denoted by v(x). Consider the problem 

' d t u(t) + A(u(t))u(t) = F(u(t)) in ft, 

Bj(u(t))=0 on 9ft, je{l,---,m}, (3.1) 

w(0) = uq in ft. 

Here we employ the maps 

[A(u)v](x) = a <*( x > u ( x )> Vm(x), • • • , V 2 ™" 1 ^)) D a v{x), 

\a\—2m 

[F(u)](x)=f(x,u(x),Vu(x) r -' ,V im - 1 u(x)), .xeft, 
[Bj(u)](x) =bj{x,u{x),Vu(x),--- ,V m 'u{x)), x G 9ft, 

for functions u E BC 2 " 1 - 1 {Ti;C N ), and v E Wg m (fl;C N ). The numbers m, are 
integers strictly smaller than 2m, and with E — C , the coefhcients are subject 
to the following regularity assumptions 

(R) a a G C l (E x E N x ••• x E 1 * 2 '*' 1 ; BC(U; B{E))) for \a\ = 2m, 
f E C l (E x E N x •■■ x E N * m ~\BC(p,;E)), 

b , e C 2m+l- mj (QQ x E x E N x . . . x E N m i . £ ) for j g . . . 5 m }_ 

We set 5 = (B 1} ■■■ , S m ). We point out that, for a fixed u E BC 2 " 1 " 1 ^, C N ), 
A(uq) is a linear differential operator of order 2m with bounded coefficients; 
whereas F contains all terms involving derivatives of order \a\ < 2m. 

We will employ the L p -setting for this problem as in [19], hence we fix p > n+2m 
and the basic spaces 

X =L p (n;E), X 1 = W* m (Sl;E), X 7 = (X ,X 1 ) 1 _ 1/p = W^ 1 - 1 / p \ft;E). 

As in Section 2 we denote the norm in Xj by | • \j and open balls in Xj by 
Bxj(u,r), j = 0,1,7. Note that by the Sobolev embedding theorem we have 
X y > £?C 2m_1 (ft; E), which allows us to plug in functions u E A 7 into the 
coefficients of A, into / and into the functions bj pointwise, without any growth 
restrictions on these nonlinearities. 

Assume we have a C 1 -manifold of equilibria ^ : U — > X\ where U C M. k is an 
open neighborhood of 0, 

A(*(0)*(C) = ^(*(0) inn, {eu, 

B(*(0) = 0, on 9ft, C G U, ^ ' ' 



X E ft, 



(3.2) 
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and set u* = ^(0). Assume that the rank of \t'(0) is k and that there are no other 
equilibria near u* in Xi, i.e. E n Bx x {u*,r{) = ^(U), for some r\ > 0, where as 
in Section 2, £ denotes the set of equilibria of (|3.ip . 

The linearization of (|3.1[) at is given by the operator Ao defined as follows 

A*v = A(u*)v + (A' (u*)v)u* — F'(u*)v in Cl, 

B*v = B'(u*)v on 90, 

(3.4) 

with v G D(A*) = D(B*) = W p 2m (0; E), 
= ^*|at(b»)- 

Next we consider the property of maximal L p -regularity for the pair (j4*, -B*)> and 
in particular for the operator A . For this we only need to consider the principal 
parts of the corresponding differential operator and of the boundary operators, i.e. 

A # {x,D)= a Q (a: ) u,(x),...,V 2m - 1 u,(x))D a ) 

\a\=2m 

B^(x,D) = i^\db j Jd{d^u)\{x,...^-^{x))Df i , 

\(]\=m 3 

for j = 1, . . . , m. Note that we use the notation D = — iV, hence V' 3 = i^^D 13 . It is 
shown in |12] that normal ellipticity of A# and the Lopatinskii-Shapiro condition 
for (A#, B#) are necessary, and in [llj that they are also sufficient for L p -maximal 
regularity of Aq. These conditions read as follows. 

(E) For all xeQ,£e R n , |£| = 1, a(A#{x,£)) C C+, 
i.e. A(x,D) is normally elliptic. 

(LS) For all x G dQ, £ G R n , with £ ■ v[x) = 0, A G C+, A ^ 0, and h G E m , 
the system of ordinary differential equations on the half-line 

Xv(y) + A#(x, £ + iv(x)d y )v(y) = 0, y > 0, 

Bj#(x,£ + iv(x)d v )v(0) = hj, j = l,...,m, 

admits a unique solution v G Cq(M.+ ; E). 
(Lopatinskii-Shapiro condition.) 

Now assume that u* G Xi is an equilibrium of (|3.ip . and let conditions (R), (E), 
and (LS) be satisfied. It was shown in [T!5] that (|3.ip then admits a local strong 
solution in the L p -sense for each initial value uq G X 7 , provided the compatibility 
condition B(uq) — holds and \uq — it*| 7 is sufficiently small. The solution map 
[uq h- » u(t, uq)] defines a local semi-flow in X 7 near on the nonlinear phase- 
manifold 

M = {ueX 1 : B{u) = on dVl\. 

In case the equilibrium is hyperbolic, i.e. o-(A )rMM. = 0, it was moreover shown 
in [19] that it is isolated and that it has the so-called saddle point property, which 
means that the local semi-flow in M. admits a unique stable and unstable manifold 
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near u* . We refer to [19] for details as well as to [24] in the case of linear boundary 
conditions. 

Returning to our situation, differentiating (|3 . 3|) w.r.t. £ we obtain for £ — 

A(u*)*'(0) + [A'(u*)V'(Q)]u* - F'(u*)*'(0) = in Q, 

Bj(u*)*'(0) = on dft, j = l,...,m. 

This shows that the image of \&'(0) is contained in the kernel N(Ao) of Aq, and 
also that T Ut (£), the tangential space of £ at u*, is contained in 7V(Ao). As in 
Section 2 we assume now that i?(\&'(0)) = N(Aq), that the eigenvalue of Aq 
is semi-simple, and that the remaining spectrum of Aq is contained in the open 
right half- plane C+. Note that by boundedness of f2 and compact embedding, the 
spectrum of Aq consists only of isolated eigenvalues of finite algebraic multiplicity, 
anyway. We can now state the main result of this section. 

Theorem 3.1. Let 2m + n < p < oo, let D, C ffi™ be an open bounded domain 
with boundary of class C 2m , and let the spaces Xj, j = 0, 1,7, be defined as above. 
Suppose u* £ Xi is an equilibrium of H3.1\) , and assume that conditions (R), (E), 
and (LS) are satisfied. Let Aq defined in (|3.4|) denote the linearization of VS. 1\) at 
u*, and suppose that it* is normally stable, i.e. assume that 

(i) near u„ the set of equilibria £ is a C 1 -manifold in X\ of dimension k G N, 

(ii) the tanqent space for £ at it* is given by N{Aq), 

(hi) is a semi-simple eigenvalue of Aq, i.e. R(Aq) © N(Aq) = Xq, 

(iv) a(Ao) \ {0} C C+ = {z e C : Re z > 0}. 
Then u* is stable in X-y, and there exists 5 > such that the unique solution u(t) 
of \3.1}) with initial value uq £ X-y, satisfying \uq — u*| 7 < S and the compatibility 
condition B{uq) — on dfl, exists on R+ and converges exponentially fast in X 1 
to some Uoo £ £ as t — > oo. 

Proof, (a) The proof is similar to that of Theorem 12. II It is based again on the 
reduction to normal form. We use the notation introduced above and denote as 
in Section 2 by P s and P c the projections onto Xq = R(Aq) resp. X c — N(Aq). 
We first center (|3.ip around u„ by setting u — u — u*, and obtain the following 
problem for u. 

dtu + A*u = G{u) in il, 

B*u = H(u) on dD,, (3.5) 

u(0) — Uq : = Uq — u» in O. 

Here G is defined as in Section 2, and 

H(u) — B*u — B{u* + u) = — [B(u* + u) — S(u») — B'(u*)u]. 

Exactly as in the proof of Theorem l 2 . 1 1 we obtain a function <f> £ C 1 (i?x c (0, po), X[) 
with 0(0) = 0'(O) = such that the equilibrium equation 

A*V(C) = G(^(0) in n, C € U, 

b^(c) = H(i>{0) on an, c e u, 
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for (13. 5|) can equivalently be expressed by 
P c A^(v) = P c G(v + cf>(v)), 

P s A^(v) = P s G(v + cj)(v)), B^(v) = H(v + (/>(v)), '"' > ' >>) 

for every v € Bx"(0, po). We can now introduce the normal form of (|3.1|) for the 
variables 

v := P c {u - u„) = P c u, w := P s (u - u„) - 0(P c (u - u*)) = P s u - 4>{P c u), 
which reads as 

d t v = T(v,w) in f2, 

d t w + P s A.,P s w = R{v, w) in Q, 

B*w = S{v,w) ondQ, ^'^ 

v(0) — wo, w(0) = wo in f2. 

Using the equilibrium equations in (|3.6[) we can derive, similarly as in Section 2, 
the following expressions for T, P and S: 

T(v, w) = P c (G{v + <f>{v) +w)- G(v + tf>(v))) - P c A*w, 

R(v,w) = P s (G(v + (j)(v) + w) - G(v + cf)(v))) -(f>'(v)T(v,w), 

S(v, w) = H(v + (j)(v) +w)- B*<f){v) 

= H(v + (j)(v) + to) - H(v + tf>(v)). 

Clearly, 

R(v, 0) = T(v, 0) = S(v, 0) = 0, ve B xc (0, p ). 
Therefore we are in the same situation as in Section 2, except that here the infinite 
dimensional part, i.e. the equation for to, has a nonlinear boundary condition in 
case S(v, w) ^ 0. 

(b) Let < a < oo and define the following function spaces on (0, a) x 57: 

Ei (a) = H^((Q,a);Xa)nL p {(Q,a);X 1 ), E (a) = L p ((0, a);X ). 

We also need spaces for the boundary values. For this purpose, we set with «y = 
1 — rrij/2m — l/2mp 

Y = L p (dn ; C N ), Yj = wl mK > (dn ; C N ), 

and 

m 

F(a) = JjF^o), W j (a) = W^((0,a);Y o )nL p ((0,a);Y j ), j = l,...,m. 

3 = 1 

Note that by trace theory we have 

¥ 3 (a) ^ BC([0, a}; W^ nK ^ 2m/p (dfl; C N ), 

W^ mKJ - 2m/p (dil;C N )) ^ BC(dQ;C N ) 

by the condition p > 2m + n and since nij < 2m. The spaces Fj-(a) are the trace 
spaces on the lateral boundary (0, a) x dfl of (0, a) x for the derivatives D^u of 
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order |/3| = mj for u G Ei(et). 

The basic solvability theorem for the fully inhomogeneous linear problem 

!d t u + A*v, = f(t) in ft, t>0, 
B*u = g(t) on dfl, t > 0, (3.8) 
u(0) — uq in 
in the L p -setting reads as follows, see [12] . 

Proposition 3.2. Let a < oo. The linear problem \3. 8\) admits a unique solution 
u G Ei (a) if and only if f G Eo(a), g G F(a), uq G X 7) and the compatibility 
condition B*u — g(0) holds. There is a constant C = C(a) > such that the 
estimate 

Mkfa) < C(\U \ 7 + ||/||kp (a) + ll3llF(a)) 

holds for the solution u of \3.8\) . 

We shall also need a variant of Proposition 13 . 21 for the problem 
' d t w + P s A*P s w = f{t) in ft, t > 0, 

B*w = g(t) on dfl, t > 0, (3.9) 
w(0) — wo in fl 

on the half- line, where we assume wo G and / G L P (M.+ ; Xq). For this purpose 
we proceed as follows. Suppose first that u solves (13.8[) with uq — wq. Since 
A*P c u = B*P c u = we then conclude that w — P s u solves problem (|3.9p . Let U\ 
denote the solution of (|3.8[) with A* replaced by ^4* + 1 . The spectrum of A + 1 is 
contained in C+, hence we may apply Proposition 3.1 of [19) to obtain a uniform 
estimate for u\ in Ei(oo). Then U2 = u — u\ solves the problem 

d t u 2 + A*u 2 = mi, B*u 2 = 0, tt 2 (0) = 0. 

As o~(A s ) C C+, A s has maximal L p -regularity on the half-line, hence we obtain 
also a uniform estimate for P s u 2 in Ei(oo). These arguments yield the following 
result. 

Proposition 3.3. Let a < oo. The linear problem \3. 9\) admits a unique solution 
w G Ei (a) n L p ((0,a);X$) if and only if f G L p ((0, a);X(j), g G F(a), w G X*, 
and the compatibility condition B*wq — g(0) holds. There is a constant Cq > 0, 
independent of a, such that the estimate 

Hk(a) < C (|w | 7 + |/|k(a) + \\gh(a)) 

holds for the solution w of 13. 9\) . for all functions f G L p ((0, a); Xq), g G F(a) and 
all initial values wq G X^. 

Proposition 13.31 remains valid when we replace w(t) by e at w(t), f(t) by e at f(t) 
and g(t) by e at g(t) where < a < uj, uj < inf{ReA : A G cr(A s )}. 

(c) Next we consider the nonlinearities R, T, and S. Since by assumption the 
functions a a and / are in C 1 and p > n + 2m, it follows easily via the embedding 
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X-y <^-> BC(Vt) that A and F are as in Section 2. Hence we obtain as there the 
estimates 

|T(w,u/)|o<C 2 (^ + r)Hi + C'3Mi, «eB Xc (0,p), u>e5 x ,(0,p)nX l7 

where r = 3p and C3 := |P c >l.*P c |Bnti,;c<=\. Since <fi'(0) = we can assume that 
po was chosen so small that W [ w )\b{X" ,x s ) < V f° r au w G Px°(0,/Oo). With this 
we obtain 

\R(v,w)\ < dfa + rJMi, « G B X o(0,p), w G B x «(0,p)nJCi. 

Observe that in contrast to the previous section the constant C3 is no longer 
small since P l , I G {c, s}, and ^4* do not commute. However, this does not alter 
our conclusions. It is more involved to derive the estimates on S needed for 
Proposition 13.31 Fortunately, we can refer to [HI Proposition 3.3]. This result 
implies 

|e wt (ff(«i) - fl"(ua))|F(a) < V\\e ut (ui ~ « 2 )k(a), 

for all e ut ui,e ut v,2 G Ei(a) such that |ui(i)| 7 , |u 2 (i)| 7 < r, t G [0,a]. Therefore, 
by possibly decreasing r > 0, with U\ — v + 4>{v) + w and u 2 — v + <fi(v) this yields 

\\e ut S(v,w)\\ ¥ ( a) < ^le^iolE^o), 

for all e^v^^w G E x (a) such that u([0, a]) C B X c(0,p) andw([0,a]) C B x -(0,p). 
These are the estimates we need for applying Proposition 13.31 

(d) We may now follow part (c)-(f) of the proof of Theorem 12. II to complete the 
proof of Theorem 13. 1[ the needed local well-posedness result being Proposition 4.1 
in [19]. □ 

4. Convergence of solutions for the Mullins-Sekerka problem 

We consider the two-phase quasi-stationary Stefan problem with surface ten- 
sion, which has also been termed Mullins-Sekerka model (or Hele-Shaw model with 
surface tension) and is a model for phase transitions in liquid-solid systems. Let 
Q be a bounded domain in K™, n > 2, with smooth boundary dil. Let F C ft 
be a compact connected hypersurface in f2 which is the boundary of an open set 
fio C ft, and let T{t) be its position at time t > 0. Denote by V(t, •) and n(t, ■) the 
normal velocity and the mean curvature of T(t), and let £ti(t) (liquid phase) and 
f2 2 (t) (solid phase) be the two regions in SI separated by T(t), with Sli(i) being 
the interior region. Let further v(t, •) be the outer unit normal field on T(t) with 
respect to fii(i). We shall use the convention that V is positive if fii(t) is expand- 
ing, and that the mean curvature is positive for uniformly convex hypersurfaces. 
The two-phase Mullins-Sekerka problem consists in finding a family T(t), t > 0, of 
hypersurfaces satisfying 



V=[d v u K ), 00, r(0) = F o , 



(4.1) 
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where u K = u K (t, •) is, for each t > 0, the solution of the elliptic boundary value 
problem 

Au = in fi x (t) U S7 2 (t) 5 

u = k onr(t), (4.2) 
<9„u =0 on dVl. 

Here [c^w^] := d v u 2 K — d v u\ stands for the jump of the normal derivative of u K 
across the interface r(i), and d v u denotes the normal derivative of u on dCl. 

Assuming connected phases and that the interface does not touch the fixed 
boundary dfl, the set of equilibrium states of (|4.ip . (|4.2p consists precisely of all 
spheres Sr(xq) C ft, where R denotes the radius and Xq the center. Thus there is 
an (n + l)-parameter family of equilibria, the parameters being the n coordinates 
of the center xq and the radius R. 

Let now S C be some fixed sphere without boundary contact. We are in- 
terested in the asymptotic properties of solutions of the Mullins-Sekerka problem 
that start in a neighbourhood of E, that is Tq is close to E. Following [15] we first 
use Hanzawa's method to transform the original problem to a system of equations 
on a fixed domain. Here the basic idea is to represent the moving interface T(t) 
as the graph of a function in normal direction of a fixed reference surface, which 
will be S in our case. Denoting the parameterizing function by p(t, ■) this leads to 
a problem on E of the form 

p + B(p)S(p) = 0, t>0, p(0) = p , (4.3) 

where S(p) is the solution of the transformed elliptic boundary value problem 

' A(p)v =0 in fix U Q 2 , 

v = K(p) on £, (4.4) 

d v v =0 on 9f2. 

Here Oi and are the two regions in separated by S, with f^i being enclosed by 
E. By construction, the solution T(-) = E of the original problem (|4.1|) corresponds 
to the solution p = of ([43]) . The operator F(-) := B(-)S(-) in (|33]l is a nonlocal 
pseudo-differential operator of third order and renders (|4.3|) a quasilinear parabolic 
problem, see [15] for its precise definition and more details. 

We want to study (|4.3[) in an L p setting. Let p > n + 2 and define 

X = X 1 = W*- l 'v (E). 

Given J — (0,a), a > 0, we view (|4.3p as an evolution equation in the space 
Eo(J) = L p (J;Xq), that is we are interested in solutions of (|4.3[) in the class 
Ei (J) = Hp(J; X ) n L p (J; X{). For the corresponding trace space we have 

X, := (Jfo.^Oi-i/pj, = <- 4 /p(E). 

Note that, by Sobolev embedding, we have X\ <—* X 7 C 2 (E). 
Let 

21 := {p e C 2 (E) : |p| c(E) < V } 
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with 77 > sufficiently small denote the set of admissible parameterizations. Set- 
ting U := Xi n 21, one has F € C°°(U;X ) and the linearization L := F'(0) is 
given by 

Lp = -[dvTAvp], 
where Tg denotes the solution of the elliptic problem 

Av = in fii U fi 2 , 

u = .9 on E, (4.5) 
d v v = on 90, 
and 

1 m-i 

Ay, : 



n - 1 V i? 2 

with i? being the radius of E and the Laplace-Beltrami operator on E. The 
operator Ay is the linearization K'(0) of the transformed mean curvature operator 
K(p) at p = 0. Concerning the linearization L = F'(0) we refer to [THIUS]- 

One can show (cf. [551 Theorem 2.1]) that the spectrum of L consists of count- 
ably many real nonnegative eigenvalues of finite algebraic multiplicity, and that 
is a semi-simple eigenvalue of L with multiplicity n + 1, see also [TSJ Propo- 
sition 5.4 and Lemma 6.1]. Moreover, the kernel of L is given by N(L) = 
sp&n{Y ,Yi, . . . ,Y n }, where Yo = 1 ; an d where Yj, 1 < j < n, are the spheri- 
cal harmonics of degree 1. We may assume that Yj = R~~ Pj\s> 1 < j < n , with 
Pj being the harmonic polynomial of degree 1 given by pj (x) — Xj for x € K.™ ; by 
Pj\y we mean the restriction of pj to E. 

Let us assume that E is centered at the origin of R n . Suppose S C fl is a sphere 
that is sufficiently close to E. Denote by {z\, . . . , z n ) the coordinates of its center 
and let zq be such that R + z corresponds to its radius. Then, by [TH Section 6], 
the sphere S can be parameterized over E by the distance function 



P (z) = J2^-R- 

3 = 1 



s vZ ^) 2 + (*+*b) a -2.4 



Denoting by O a sufficiently small neighbourhood of in W l+ , the mapping 
[z h- > p(z)] : O — > Vt^p 1 ^ P (E) is smooth and the derivative at is given by 



p'(0)h — hjYj, heR n+1 . 
3=0 

So we see that near E the set £ of equilibria of (|4. 3[) is a smooth manifold in Xi 
of dimension n + 1, and that the tangent space Ts(£) coincides with N(L). 

In order to be able to apply Theorem 12.11 from Section 2 it remains to verify 
that the operator L has the property of maximal L p -regularity. This means we 
have to show that for any J = (0, a), a > 0, and any g £ Eo(J) the problem 

p + Lp = g, teJ, p(0) = 0, (4.6) 
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has a unique solution in the space Ei(J). By means of the standard localization 
method, perturbation arguments, and by solving certain elliptic auxiliary prob- 
lems, (|4.6I) can be reduced to the following two-phase problem on R™ x R with 
R = R\ {0}: 



~A x w 




= 0, 


te J, x e R", y e R, 


(4.7) 




\-A x a 


= 0, 


teJ,xe R", 


(4.8) 


d t <J - 


[dyW] 


= h, 


teJ,xe R", 


(4.9) 




<r(0) 


= 0, 


sc 6 R n . 





Here [9 y w] = 9 y w| y=0 + — 9 y u;L = o-, an d ft S Lj,(J; Wp _1//p (R n )) is a given func- 
tion. We take the Fourier transform w.r.t. x and denote the transformed functions 
by w and a. Then (|4. T[) and (I4.8[) imply that w = e - ^' \£\ 2 <T- Inserting this into 
(|4.9[) leads to the subsequent problem for a on R": 

d t d + 2\^a = h, te J, cr(0) = 0. 

Set r = Wp 1_1/p (R") and Yi = Wp 4_1/p (R n ) and let G be defined by G = d/di 
with domain D(G) — QHp(J;Yo), here the zero means vanishing trace at t = 0. 
Then G is sectorial, invertible and admits an 7^°°-calculus in L p (J;Yq) of angle 
ir/2. Let further D be the operator in L p (J;Yq) with symbol 2|£| 3 and domain 
D(D) = L p (J;Yi). Then D is sectorial and admits an 7i°°-calculus in L p (J;Yq) 
of angle 0. Thus by the Dore-Venni theorem, the equation Gtr + Da = h possesses 
a unique solution a £ -D(G) n D(D). Hence L has the property of maximal L p - 
regularity. 

So all assumptions of Theorem 12.11 are satisfied, hence we obtain the following 
result, which is the main result in [15] except for the different functional analytic 
setting. 

Theorem 4.1. Let p > n + 2 and £1 C R n be a bounded domain with boundary 
of class G 2 . Suppose E is an arbitrary sphere in fl of radius R without boundary 
contact. Then p = is a stable equilibrium of in X~ t — Wp 4 ^ P (S), and there 
exists 5 > such that if \po\-y < 5, then the corresponding solution of exists 
globally and converges at an exponential rate in X y to some equilibrium as 
t — * oo. In this sense, the sphere S is a stable equilibrium of the Mullins-Sekerka 
problem, and any solution T(-) of l{4.1\ ) that starts sufficiently close to E exists 
globally and converges to some sphere at an exponential rate as t — > oo. 

This approach can also be used to show the stability of spheres for the two-phase 
quasi-stationary Stokes flow in a bounded domain, see [17\ 116] for alternate ap- 
proaches in the one-phase case. Moreover, it can be applied to models in tumor 
growth, see [3] for a discussion of existing work. 
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5. Stability of travelling wave solutions to a quasilinear parabolic 

equation 

The situation of the generalized principle of linearized stability may occur when 
studying the stability of travelling wave solutions of parabolic equations, see e.g. 
[18] Section 5.4] for the semilinear case. In what follows we want to consider a 
quasilinear variant of the Huxley equation: 

u t - (£r(«x)) a = /(«). * > 0, i£E. (5.1) 

Here f(r) = r(l — r)(r — a), r 6 R, where a G (0, 1/2) is a constant, and a is a C 2 
smooth function on R satisfying 

< ci < cr'(r) < c 2 , re R. (5.2) 

A travelling wave u(t, x) = iu(x + Vt) with speed V satisfies 

(ct(w'(s)))' -Vw'(s) + f(w(s)) =0, set. (5.3) 

Similarly to the special case a(r) = r (cf. [IE]), one can show, by means of a phase 
plane analysis, that for some V > (|5.3|) admits a solution w with — * as 
s —> — oo and w(s) ~> 1 as s — * oo. For this purpose we introduce the variable 
z := w'. Then (|5.3p is equivalent to the system 

; = z, 

■>. = ^(v z -fM). 

Denoting its right-hand side by H(w, z), we find that 



H'(i,0) 



1 

f'(i) V 



T'(0) CT '(0) 



it {0,1}. 



We have /'(0) = —a and /'(l) = a — 1, thus the equilibria (0,0) and (1,0) are 
both saddle points. The eigenvalues of H'(0,0) are given by 

(1/Ai,l) is an eigenvector to Ai > 0, thus the unstable manifold points into the 
first and third quadrant, with steeper slopes for higher values of V > 0. 
Define the functions 

F(V) = f V f(r)dr, G(y) = f a'(r)r dr, y e R. 
Jo Jo 



Then (|5.3p implies that 
d 



— (g(z(s)) + F{w{s))) =Vz(s) 2 , set. (5.4) 

In particular G(z) + F(w) is a first integral if V — 0. By (|5.2p we further have 
Ciy 2 < 2G(y) < C7 22 / 2 , y e M. 
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We now consider the trajectory 7 that (near the origin) lies on the unstable 
manifold to (0, 0) in the first quadrant. For V — 0, 7 cannot reach the line w = 1, 
since F(w) < on 7, and F(l) = ±(± - a) > 0. In case V > 0, (|5.4|) shows that 
7, as long as it remains in the first quadrant, moves through increasing values c of 
the level curves G(z) + F(w) = c. For V sufficiently large, 7 will reach the level 
curve to c = F(l) at some point with z > 0. For continuity reasons, there exists 
then V > for which 7 becomes a heteroclinic orbit, connecting (0,0) and (1,0); 
observe that 7 is the only such orbit. Hence there is a smooth solution w to (|5.3p 
satisfying w'(s) > for all s £ R and 

(w(s), u>'(s)) — ► (0,0) as s -4 —00, (w(s),w'(s)) — > (1,0) as s — > 00, (5.5) 

both exponentially fast. Clearly each translate w(- + a) with a £ R enjoys the 
same properties. 

In order to investigate the stability of the travelling wave we change to moving 
coordinates with y := x + Vt and u(t, y) :— u(t, x), thereby transforming (|5.ip into 

u t - (a(u y ))y + Vu v = f(u), t > 0, y £ R. (5.6) 

Evidently, u = w and all translates of it are equilibria of (|5.6|) . Letting v :— u — w 
be the deviation from w, the equation for v reads as 

v t - (<r(v v + w y ))y + V(v y + w y ) = f (v + w) , t>0, J/ g f . (5.7) 

For this equation, the set 

£ = {w(- +a) - w(-) : a £ R} (5.8) 

forms a one-dimensional smooth manifold of equilibria. Observe that £ contains 
all equilibria w(y) of (|5.7[) satisfying w(y) — > as \y\ — > cxo. By definition, the 
travelling wave under study is stable, if v = is stable for (|5.7[) . Of course, this 
has to be understood in the sense of a suitable functional analytic setting. 

Let us choose again the L p setting to study (|5.7|) . Let 1 < p < oo, X = L p (R), 
and Xx = H%(R). Define M : X x -> X by 

M(«) = -(cr(w y + w a )) a + t/(w y + w a ) - /(« + w), V £ X 1; 

This definition makes sense since by (I5.2I) - (|5.3[) and (|5.5p w^y, and /(u>) belong 
to L P (R). The fact that f(v + w) belongs to the space L p (R) can be justified by 
observing that 

f(v(y)+w(y)) = f(w(y))+[ f'(w(y) + rv(y)) dr v(y), y £ R. 

Jo 

Since v and w belong to Cq(R) one readily verifies that f'(v + w) is continuous and 
bounded, and this yields the statement for f(v+w). The linearization Ao :— M'(0) 
of M at v = is given by 

A v = -(a'(wy)v y )y + VtJj, - /'(w)w, « £ £>(A ) := X x . 

By (|5.2p , Ao is a uniformly elliptic operator with smooth coefficients whose leading 
coefficient satisfies —a'(w'(y)) — » — cr'(0) as |y| — > oo. Thus, by [TTJ Theorem 5.7], 
enjoys the property of maximal L p -regularity. 
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Next observe that £ C i?~ (R) Xl and that the tangent space for the manifold 
£ at v = coincides with spanjw/}, the span of to' G Xi. On the other hand, by 
differentiating (|5.3[) we see that >lo w ' — 0. So to show normal stability of v = 0, 
it remains to prove that is a simple eigenvalue of A and that the remainder of 
the spectrum of A lies in {z £ C : Rez > 0}. We proceed similarly as in [T51 p. 
131], generalizing the proof given there to the quasilinear case. Suppose A with 
Rc A < is an eigenvalue of Aq with eigenfunction v G X\ , that is 

^' + W + 4^, + -^ = o, yeR, 

a'(w'(y)) a'(w'(y)) 

where 

2cr'(w'(y)) 

By studying the characteristic equation for the limits y — > ±oo, one sees that there 
exist constants 8, C > such that 

\v(y)\ < CeV+^'Wy, y<0, \v(y)\ < Ce~ Sv , y > 0. 

Letting <p(y) = v(y) exp(/* b(r) dr), we have at least tp(y) = 0( e - y l^/( 2CT '(°))) for 
\y\ —>■ oo, and 

^ + ( ;? ^ ^ + - - b ^ 2 ) *> = °' f e M - ^ 5 - 9 ) 

Vo-'(to'(y)) o-'(u/(y)) / 

The function ^(y) := w'(y) exp(J Q y 6(r) dr) is strictly positive, and since Aqw' = 0, 
it satisfies 

7'Mw)) 

'K(y)) 
Combining ([BT9| . (ISTO) yields 



Vcr' w' y / 



A V" 
cr'(w'(y)) V 

Multiplying this equation by integrating over R, and integrating by parts gives 



*{y? { [(Re V {y)im)'} ' + [(Im ¥■(»)/«(»))']*} *• 



hence A = 0. In view of (|5.9p we may then suppose that ip is real and nonnegative, 
and so the last formula above implies that tp/ip — v/w' is constant on R. Hence 
N(Ao) = span{w'}. 

We now show that the essential spectrum <j e (Ao) of Aq, that is the set of all 
spectral points except isolated eigenvalues of finite multiplicity, is contained in 
[Re z > 0]. This can be seen as follows. 

By the asymptotics of to, see (|5.5p . we know that f'(w(y)) — > /'(0) = —a as 
y — > — oo, and that f'(w(y)) — » /'(l) = —(1 — a) as y — > oo at an exponential rate. 
Due to a 6 (0, 1/2), a is the smaller of the two numbers {a, (1 — a)}. Fix e > 
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so that a — e > 0. We can then find a number R > so that —f'(w(y)) > a — e 
whenever \y\ > R. Let c be a bounded continuous function on R that agrees with 
—f'(w) on [\y\ > R] and satisfies c(y) > a — e for all y 6 R. Let B : X\ — > Xo be 
the operator defined by 

B« := —{& '{Wy)Vy)y + V^Wy + C?7 . 

One readily shows that B is accretive on L P (R) and also that u(B), the spectrum of 
B, is contained in [Re z > (a—e)]. Next, note that Aq can be written as Aq — B+S, 
where S is a perturbation which is relatively compact with respect to B. For this 
we note that S can be written as Sv = —(f'(w)+c)xv, where x is a smooth cut-off 
function for the interval [-R, R], with support contained, say, in f2 = (— 2R, 2R). 
Since ffp(fi) is compactly embedded in L p (Q) for s > we conclude that S is a 
compact operator from ff5(R.) — » L p (M). A well-known perturbation result now 
shows that a e (Ao) must also be contained in [Rez > (a — e)]. Since this is true 
for every e, we have proved that a e (Ao) C [Rez > o], Aq might still have isolated 
eigenvalues outside of this set. As Aq generates an analytic semigroup, they must 
be contained in an appropriate sector with opening angle 9 < n/2. Since we have 
already shown that there are no eigenvalues in [Re z < 0] we conclude that there 
is a number a > so that <t(Aq) \ {0} C [Rez > a). Finally, since the operator 
Ax defined by 

Aiv := v yy + (^M - b \y) - b(y) 2 )v, v € D(Ai) = X u 

is self-adjoint in L2(M), it follows in view of (|5.9[) with A = that the eigenvalue 
of Aq is semi-simple. 

Summarizing we have shown that 

• the set £ consists of all equilibria of (|5.7p in X\ and forms a smooth 1- 
dimensional manifold, 

• T a (£) = span{w , } = N(Ao) 1 

• the eigenvalue is semi-simple, 

• a(Ao) \ {0} C C + ={z£C: Re z > 0}, 

and we can now formulate our main result for this section. 

Theorem 5.1. Let p G (1, oo), f(r) = r(l - r)(r — a), r € K, with a e (0, 1/2), 
and a £ C 2 (R) such that i5.ty) holds. 

Then 15.1]) possesses a travelling wave solution u(t, x) = w(x + Vt) with speed 
V > and profile w G C 3 (M) satisfying w'(r) > 0, r £ R, and (w(r),w'(r)) — » 
(0,0) as r — > — oo, as we// as (w(r),w'(r)) — > (1,0) as r — * oo. 

XTws travelling wave is stable in the sense that v = is a stable equilibrium of 
\5. 7| ) in X-y — Wp 2 ^ P (R). Moreover there exists 5 > such that if \va\-y < 5, 
then the solution v of \5. 7j ) with v(0) — vq exists globally and converges at an 
exponential rate in X y to some equilibrium Woo, i.e. to some element of the set 
£ defined in i5.8\) . In this sense, any solution u of \5.1}) that starts sufficiently 
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close to w exists globally and converges at an exponential rate as t — > oo to some 
translate w(x + Vt + a), o£l, of the travelling wave solution. 

This approach applies to many other travelling wave solutions of quasilinear par- 
abolic systems, as soon as the condition of normal stability is satisfied. 

6. Convergence for abstract quasilinear problems II 

We return to the setting of Section 2 for the case that o-(Aq) also contains an 
unstable part, i.e. we now assume that 

a(A ) = {0} Ucr s Uer u , with a s C C+, a u C C_, (6.1) 

such that o~ u ^ 0. In this situation we can prove the following result. 

Theorem 6.1. Let 1 < p < oo. Suppose u* € V D X\ is an equilibrium of 12.1]) . 
and suppose that the functions (A, F) satisfy (|2.2|) . Suppose further that A(u*) 
has the property of maximal L p -regularity. Let Aq be the linearization of \2.1]) at 
it*. Suppose that u„ is normally hyperbolic, i.e. assume that 

(i) near it* the set of equilibria £ is a C 1 -manifold in X\ of dimension m G N, 

(ii) the tangent space for £ at it* is given by N(Aq), 

(iii) is a semi-simple eigenvalue of Aq, i.e. N(Aq) © R(Aq) = Xq, 

(iv) a(A ) F\iR = {0}, a u := a(A ) n C_ ^ 0. 
Then u* is unstable in X 7 and even in Xq. 

For each sufficiently small p > there exists < 5 < p such that the unique 
solution u(t) of \2.1]) with initial value uq £ Bx^{u*,8) either satisfies 

• distx.y(w(to)> £) > P for some finite time to > 0, or 

• u(t) exists on M + and converges at an exponential rate to some Uqo 6 £ 
in X~f as t — > oo. 

Proof. The first assertion follows from [2U Theorem 6.2], so we need to prove the 
second claim. 

(a) Let P l denote the spectral projections corresponding to the spectral sets 07, 
where a c = {0} and a s ,a u are as in (|6.1j) . Let Xj = P l (Xj), I € {c,s,u}, where 
these spaces are equipped with the norms of Xj for j 6 {0, 1, 7}. We may assume 
that X\ is equipped with the graph norm of Aq, i.e. \v\i := \v\q + |Ao«|o for 
v G X\. Since the operator — Aq generates an analytic Co-semigroup on Xq, o~ u 
is a compact spectral set for Aq. This implies that P u {Xq) C X\. Consequently, 
Xq and X" coincide as vector spaces. In addition, since A u , the part of Aq in 
Xq , is invertible, we conclude that the spaces X™ carry equivalent norms. We set 
X u := X% = XI and equip X u with the norm of X , that is, X u = (X u , | ■ | ). 
As in the proof of Theorem 12.11 we obtain the decomposition 

Xi = X c © X[ © X u , Xq = X c ®X s q®X u , 

and this decomposition reduces Aq into Aq ~ A c © A s © A u , where Ai is the part 
of Aq in X Q for / e {c, s, u}. It follows that o-(Ai) = ai for / 6 {c, s, u}. Moreover, 
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due to assumption (iii), A c = 0. In the sequel, as a norm in Xj we take 

\v\j = \P c v\ + \P s v\j + \P u v\ for j = 0, 7, 1. (6.2) 

We remind that the spaces Xj have been given the norm of Xq for I e {c, u}. 

We also fix constants lo 6 (0, inf Re a(—A u )) and M 5 > such that \e Al,t \ < 
M 5 e~ ut for all t > 0. Wlog we may take w < 1. 

(b) Let $ be the mapping obtained in step (b) of the proof of Theorem 12. 1[ and 
set (j>i{x) := P l $(x) for / E {s,u} and for x £ B X a(0,p o ). Then 

<t>i €C\Bx*{Q,po),X[\ ^,(0) = $(0) = for Ze{s,u}. (6.3) 

These mappings parametrize the manifold £ of equilibria near u* via 

[x^ (x + cj) s (x) + 4> u (x) + «*)], a; e Bxc(0,p )- 

We may assume that po has been chosen small enough so that 

<1, x£Bx40, Po ), le{s,u}. (6.4) 

(c) The equilibrium equation (|2.7p now corresponds to the system 

P c G(x + (j> 8 {x) + (f> u (x)) = 0, 

(6.5) 

P l G{x + <p s {x) +<(>u{x)) = Ai^ix), xeB X c(0, Po ), le{s,u}. 
The canonical variables are 

x = P c v, y = P s v-0 s (x), z = P u v - <f> u (x) 
and the canonical form of the system is given by 

x = T(x,y,z), x(0)=x , 
y + A s y = R s (x,y,z), y(0)=y , (6.6) 
z + A u z = R u (x,y,z), z(0)=z . 

Here the functions T, R s , and R u are given by 

T(x, y, z) = P c (G(x + y + z + <j> a [x) + <p u {x)) - G(x + <f> a (x) + M%))) , 
Ri (x, y, z) = P l (G(x + y + z + 4> s {x) + <f> u (x)) - G(x + S (x) + <j) u (a;))) (6.7) 
- <t>'i(x)T(x,y,z), le{s,u}, 

where we have used the equilibrium equations (|6.5[) . Clearly, 

Ri(x,0,0) =T(x,0,0) =0, xeB X c(0 :Po ), 1 & { S:U } 1 

showing that the equilibrium set £ of (|2 . 1 1) near it* has been reduced to the set 
B X c(0 lPo ) x {0} x {0} C X c x X s x X u . 

There is a unique correspondence between the solutions of (|2.1|) close to in X-y 
and those of (|6.6p close to 0. We again call (|6.6p the normal form of (12. ip near its 
normally hyperbolic equilibrium u*. 
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(d) The estimates for Ri and T are similar to those derived in Section 2, and we 
have 

\T(x,y,z)\, \Ri(x,y,z)\ <P(\y\i + \z\), (6.8) 

for all x, z G B l x (Q,p), I G {c, u}, and y G Bx*(0,p) flXi, where p < po, r = 5p, 
and (3 = Gi (j] + r). 

(e) Let us assume for the moment that p is chosen so that 4p < pq. Let u(t) = 
+ Q(x(t)) + y(t) + z(t) be a solution of (|6.6[) on some maximal time interval 

[0,i») which satisfies distxAu{t), £) < p. Set 

ti := h(x ,y ,z ) := sup{t G (0, t«) : |u(r) - < 3p, r G [0, t]} 

and suppose that ii < i*. Assuming wlog that the embedding constant of X± <—* 
X 1 is less or equal to one it follows from ()6.2j) . (16. 4|) and the definition of ii that 

|z(t)|,|y(i)| 7 ,Kt)| <3p, tG[0,ti], (6.9) 

so that the estimate (|6.8p holds for (x(t), y(t), z(t)), t G [0,ti]. 

Since £ is a finite-dimensional manifold, for each u G -Bjj^ (u*, 3p) there is u G £ 
such that disty (u, £ ) = jit — zt| 7 , and by the triangle inequality u G Bx^ (w*, 4p). 
Thus we may write it = u* + + y + z and u = u* + $(5), and therefore 

p > distx 7 (u, 8 ) = |u — «| 7 

= |X — x| + | J/ + a (x) - 0a(x)| 7 + |z + u (x) - 4>u(x)\ 

>\x-x\ + \z\ - \<j) u (x) - (f> u (x)\ > \z\, 

since x,x G -Bx c (0,po) and 4> s is non-expansive, see (16. 4p . Therefore we obtain 
the improved estimate \z(t)\ < p for all t G [0, t{\. 

We begin the estimates with that for the unstable component z(t). Integrating 
the equation for z backwards yields 

rti 

z (t) = gfa)- J e A ^ s -^R u {x(s) iy (s),z{s))ds. (6.10) 

With dnH} and \z(ti)\ < p we get 

\z(t)\ < M 5 e-»^-Vp + (3M 5 J* e — (—*) ( | y(s) | 1 + | z(s) |) ds 

for t G [0,<i]. Gronwall's inequality yields 

\z(t)\ < Mse-^-Qp + (3M 5 J 1 e -« ^(s)^ d s 

for t G [0,ii], where u)\ = u) — /3-Ms > provided (3, i.e. 77, r are small enough. In 
particular, with Me = M^/coi, this inequality implies 

\zh t (.Ji;X ) < M eP + PM 6 \y\ L9 (_j ViXl ), (6.11) 

where we have set J\ = (0,fi); here q G [1, oo] is arbitrary at the moment. A 
similar estimate holds for the time-derivative of z, namely 

l&KfriXu) < (IKI + 0)\\4L q (JuX o ) + Pbh^Jr^y (6.12) 
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Note that 

\z{t + h)-z{t)\<h^'\z\ Lp{Jl , Xo) , 
r ti-h (6.13) 

\z(t + h)-z(t)\dt<h\z\\ LliJl . Xo y 



Next we consider the equation for x. We have 

\x(t)\<\x \+ f \x(s)\ds=\x \+ [ \T{x{s),y(s),z(s))\ds 



< \xo\ + P{\v\l 1 (j 1 ,x 1 ) + \\4l 1 (J 1 ;X ))- 

Combining this estimate with that for z we obtain 

sup \x(t)\ < \x \ + \\x\\ Ll{Jl .,x ), 
te.Ji 

\x\l,{ Jv ,x ) < P{M 6 p + (1 + f3M 6 )\\y\\ LqiJl . Xl) ). 

This estimate is best possible and shows that in order to control \x(t)\ we must be 
able to control \\y\\L 1 (j 1 -,x 1 )- Note that 

\x{t + h)-x{t)\<h l IP'\x\ Lp{JvtXo) , 

-h (6.14) 
\x(t + h)- x(t)\dt < h\\±l Ll ( Ji;Xo y 

Now we turn to the equation for y, the stable but infinite dimensional part of the 
problem. As in the proof of Theorem 12.11 part (e) , we obtain from (JOT 



IMk( tl ) < Mi\y h + PM (\\y\\ El{tl) + Mk(ti))- 

Employing (|6.11|1 with q = p we get 

IMk( tl ) < M^yo^ + (3M M 6 p + [3M (1 + pM 6 )\yl Mtl) . 

Assuming (3M Q (1 + /3M 6 )) < 1/2, this yields 

IMk(tO < 2Mi|y | 7 + 2(3M M 6P . (6.15) 

Repeating the estimates leading up to (12.24|) with a = we now get 

\y(t)\ 7 <C 5 (\y \ 7 + Pp), te[o,h], (6.16) 

where C5 is a constant independent of /?, yo and t\. In particular, we see that 
I2/WI7 < P for all t s Ji, provided \yo\-y and (3, i.e. 77 and r are sufficiently small. 



For later purposes we need an estimate for \y(t + h) — y(t)\ 7 . We have 



\y(t + h)- y(t)| 7 < C\y(t + h) - lo~' r + h)- y(t)\J 



< ch^v'iyi^j^m + + (6 ' 1?) 

for alH G [0,h], t+fc 6 [Q,*i] with y(t+h),y(t) 6I1. We remind that 7 = 

Unfortunately, this is not enough to keep \x(t)\ small on J\, for this we need to 
control |y|_Li(Ji;Ai)i and we cannot expect maximal regularity in L\. 
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To handle we are forced to use another type of maximal regularity, 

namely that for the vector-valued Besov spaces Bf QO ( Ji; X), where a 6 (0, 1); cf. 
[23l Theorem 7.5]. Before stating the result we remind that 

\\9\\b^(J-X) ■= \g\Li(J;X) + [g]j:a,X, 

ra—h 

[g]j-, a ,x ■= sup h~ a \g(t + h) - g(t)\x dt 

0</i<min(l,a) Jo 

defines a norm for g 6 Bf^J; X), where J = (0, a). The maximal regularity 
result, which is valid for all exponentially stable analytic Co-semigroups, reads as 
follows: there is a constant M7 depending only on A s and on a 6 (0, 1) such that 
the solution y of 

y + A s y = f, teJ, y(Q)=y , (6.18) 

satisfies the estimate 

WvlBf^xi) < M 7 (\yo\ DAs ( a ,oc) + Ifh^JiXg)) • 

Note that this estimate is in particular independent of J = (0, a), by exponential 
stability of e~ Ast . Further we have j/o G l 7 flX s = Da s (l — l/p,p) DA s {a, 00), 
provided a < 1 — 1/p. Another parabolic estimate valid for (|6.18p that we shall 
make use of reads as 

Wyhf^xi) < M 8 (\yo\ DAe ( a ,oo) + |/lii(Ji;x;)), 
provided a < 1. Here the constant M§ is also independent of J = (0, a). 

We set Rx{t) = -<j)' s (x(t))T(x(t),y(t), z(t)) and recall that \^ s (x(t))\ B(x ^ Xl ) < 1 
for t £ [0, tx\. Employing the Li-estimate for z, see (|6 . 1 1 1) . yields 



l-RiliiCJx;^) < / \T(x(s),y{s),z(s))\ds < f3(\\y\\ Ll{ j uXl ) + Nl^Xq)) 
Jo 

< 0M 6 p + 0(1 + (3M 6 )\\y\\ Ll{Jl . Xl) . 
Therefore, for the solution y\ of (|6.18|) with / = R\ we obtain 

lwi|uf 00 (Ji;x 1 ) < M 8 (|y | 7 +/3M 6 p + 0(l + M 6 0) || 
Next let i? 2 (t) = P s (G($(a;) + y + z) - G($(x))). Then by estimate (|2T9|) 

I-R2|Ui(Ji;X ) < ^(ll2/||ii(Ji;Xi) + II 2 II Li (Ji ;X ) ) 

< 0M 6 p + 0(l + M 6 p)\\y\\ Ll(Jl , Xl) , 

and with some constant C6 

|i? 2 (i) - ifc(*)|o < C 6 j9(|v(t) - y(t)|i + \z{t) - z(t)\ + \x(t) - x(t)\) 

+ C 6 |y(t)|i(|y(t) - y(t)| 7 + |ac(t) - x(t)\ + \z(t) - z(t)\). 
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Hence we obtain the following estimate 

rti—h 

[ifelco < C 6 f3{[y} aA + [z] a . + [x] a . } + C 6 sup h~ a / |j/(t)|i 

0<h<h± JO 

■ {\y(t + h) -y(i)| 7 + \x(t + h) - x(t)\ + \z{t + h) - z(t)\) dt 

where we set h\ := min(l,ii) and [-] a ,j '■= Hj i; q,Xj for j = 0,1. (|6.11[) (|6. 13[) 
yields for each a G (0,1) 

[z]a,o < \\zht(Jf,x ) < C 7 {p + j3\y\ Ll{Jl . Xl )), 

with some uniform constant C7. In the same way we may estimate [a;]a,o- Next 
we have again by (j6.11|) - (j6.13ll 

sup h- a "|yW|i|z(t + / l )-z(i)|^</i 1/p '^||i|| ip ( Ji; x )||2/||L 1 (,7 1 ;X 1 ) 

0<h<h 1 JO 

< C 8 (|j/ol 7 + P)IMIza(Ji;*i) 
provided a < 1 — 1 /p, and similarly for the corresponding integral containing the 
x-difference. Last but not least, for a < (1 — 7)(1 — 1/p) we have by (|6.17|) 

sup hr a [ 1 |»(t)|i|v(* + /i)-y(t)| 7 <ft 

0</i<7ii JO 

<2C/ l ( 1 -^^" Q ||y||i7 Ji;Xo) || y || ip(Ji;Xl) || 2/ ||l i(Ji;Xi) 
<C 9 (| yQ | 7 + /3p) 2 ^||y|X l(Ji;Xl) 
< Cio((||ft)| 7 + /3p) 2 + (|w)|-r + 0p)hh 1 (j 1 ;x 1 )), 
where we used Young's inequality in the last line. 

Collecting now all terms and choosing a = (1 — j)/p' = 1/pp' , we find a uniform 
constant C\\ such that for |j/o| 7 < 5 

bbf^CJiiXf) < Cn(\yoU + + (/3 + p + 5)\\y\\ B?oa{Jl .xf)) , 

hence 

bbx^Xf) < Mb^jux*) < 2Cii(|y | 7 + /3p), (6.19) 
provided Cn(/3 + p + (5) < 1/2. Choosing now first /3, i.e. 77 and r small enough, 
and then p and S > 0, we see that |w(ii) — it*| 7 < 3p, a contradiction to £i < £*. 
As in (e) of the proof of Theorem 2.1 we may then conclude that — 00, which 
means that the solution exists globally and stays in the ball Bx^(u*,3p). 

(f) To prove convergence, let (x(t), y(t), z(t)) be a global solution of (|6.6p that 
satisfies 

\x(t)\, |y(i)| 7 , \z(t)\ < 3p, for all t > 0, 
see (j6.9[) . Similarly to the proof of Theorem 12. 11 part (e), we obtain from (|6.8[) 

l^tfkcoo) < 2Mi|»b| 7 + 2/3Af ||e" t z| Eo(oo) , (6.20) 
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where lj 6 (0,inf{ReA : A € a(A a )}) is a fixed number and f3 is given in (|2.21l) . 
Repeating the estimates leading up to (I2.24p we get 

< M 2 |t/o| 7 + 2/3coMo||e^z|| Eo(oo) , t > 0. (6.21) 

From equation (|6.10l) we infer that 

/oo 
e - A -(*--)i2 u (a:(a), y(s),z(s)) ds, t > 0, (6.22) 

since |z(ii)| < p for each t\ > and e" 4 ^* 1- *) is exponentially decaying for ti — > oo. 
Using (|6.22[) and the estimate for i? u from (|6.8|) and proceeding as in the proof of 
Young's inequality for convolution integrals one shows that 

||e wt z|k(oc) < C 12 /3(le ut y\ Elioo) + \\e^ z\\ Eo{oo) ) . (6.23) 

Making (3 sufficiently small (by decreasing rj and, accordingly, r) it follows from 
(f6~2"0)) and (f6T23|) that 

|e wt y||E 1 (oo) + |e wt «|Eo(oo) < Ci3|yo|-y 

This estimate in turn, together with (|6.2ip . implies |y(t)| 7 — > and \z(t)\ —> 
exponentially fast as t — * oo. As in the proof of Theorem 12. II part (f) we get 

/>oo 

x(t) -> := x + / T(x(s),y(s),z(s))ds. 
Jo 

This yields existence of the limit 

Uoo = + Woo := u* + lim = u» + + ^(^oo) + 0«(a;oo) G f . 

i — >oo 

Similar arguments as in Section 2 yield exponential convergence of u(t) to Wqo. □ 

A result similar to Theorem 16. II is also valid in the setting of Section 3. We leave 
the details to the interested reader. 
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